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ABSTRACT 


The propagation of a simple acoustic transient in @n isovelocity 
water layer was investigated for the purpose of studying the region 
of validity of a suggested correlation between exact solutions in 
ducts and solutions for propagation in a layer. The Laplace trans- 
form method was used to obtain approximate solutions for the acoustic 
pressure, particle velocity, and particle displacement resulting 
from a step-function input in welOCrey The radiation impedance 
characteristics of a Mylar transducer were investigated to study 
their effect on the received waveforms. Computer programs were 
written to evaluate and graph the resulting waveforms. The wave- 
forms resulting from a step-function input in velocity were observed 
and compared with the predicted waveforms. With the use of a family 
of Mylar transducers, a somewhat doubtful correlation was obtained 
in the region of validity of the theoretical solutions. Possible 
causes of the disagreement between the theory and the experiment 
were investigated. It was found that the theory is believed to be 
correct, that the experiment was not accurate enough, and that 
further investigations should be conducted to study the assumed 


pressure-release characteristics of % inch polyethylene. 
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1. INTRODUCTION 


Although most of the current investigations of the propagation 
11-13 
pee ) 


$3 


of energy in waveguides have been done for electromégnetic c¢ 


3 


some have considered the propagation of pulsed acoustic energy in 


aerecuideee em The dispersion of an ultrasonic pulse propagated 


through a waveguide of circular cross-section filled with water 


has been used for the approximate determination of the frequency 


(2) 


spectrum of the original ie Jic Meseroud eb laile have shown 


that the dispersion of 2 rectangular pulse along an acoustic wave- 
guide of rectangular cross-section can be used to calibrate the wave- 
guide so that it can be used to predict the effects of simiiar 
dispersion on the transmission of pulses of other shapes that are 

not amenable to actual computation. It has been shown, by considering 
pulse compression in an acoustic waveguide, that a frequency modulated 
input pulse can be compressed into an output pulse of shorter duration 


(3) 


and of higher peak amplitude - In addition, studies of transient 


(6,7) 


propagation in layered media have been made . Many of these 
studies have used approximation techniques, such as the method of 


stationary phase, to solve the transient propagation problem. J. R. 


(8) 


Wait has presented a summary of these techniques useful for cal- 


culating the resultant propagated waveform. Using the Laplace trans- 


Ds 9 0 : ane 
(5,9, 10) have investigated transient radiation 


form technique, others 
and scattering from cylindrical bodies. 

Recently, exact solutions for the propagation of simple transients 
in waveguides have been eipieisned ae ae it has been suggested 


that there is a correlation within a certain region of validity between 


exact solutions in ducts and solutions for the propagation of 
ee yy) 
acoustic transients in isovelocity layers p 
The purpose of this research is to extend investigation of 
the region of validity of the suggested correlation and of the 


restrictions under which the solutions describing the transient 


response of the layer are valid. 
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2. THEORETICAL DEVELOPMENT 


Development of Solutions 


The response of an isovelocity layer to simple transient 
excitations will be developed using the following assumptions: 

(a) cylindrical velocity-source of radius a, 

(b) constant layer depth d, and 

(c) pressure-release upper and lower boundaries. 

The wave equation for acoustic velocity-potential in cylindrical 


coordinates, under the assumption of radially-symmetric motion, is 
ye oe 2 
See x C oc + Se = —- cs a ae O 2.1 


where r is the radial direction and z is the transverse direction. 
The convenient separation of variables 
Dicjz,t) = Ble) Ply,t) 

casts Eq. 2.1 into 


oe it & 43 - (4s) " 
Srt te or ct ot (35) Plrt) =O 


Wr \ = 
me (2 Pt [Z(2)=0 
where the radial component of the particle velocity becomes 


ur = Z(2) 32 


and the acoustic pressure becomes 
p> -e E(t) 


For pressure-release top and bottom, we have p = 0 at z = 0 


Wns nIIc 


and z = d; therefore, Z(z) =D sin —— where ae 





Thus , 
the equation for Z(z), which gives the dependence of the acoustic 


variables on the layer depth, determines the eigenfrequencies woe 
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These are the cut-off frequencies of the various radially-symmetric 
normal modes of the layer. 
Use of the Laplace transform method to solve Eq. 2.2 with 
the boundary condition at the source written as 
U, (a, ©) = 2a) V(t), 
which implies excitation of only one mode, casts Eq. 2.2 into 


the simpler differential equation 
pa — 
he ee ee 
“ee ee eee P(r, a=O ore 


where (r,s) is the Laplace transform of Mir, t),. V(t) Z(z) is 
the radial velocity at the source. The relevant boundary conditions 


become 


( 


a 


P = V(x) 


dr tr=o 
and 
gO as [ —P OO. 
The change of variable u sie | we + e Gasts Eq. 223) imto the 
form 


zo L 
——eee to 4 


u& 


= 1S 
Sl 
O 


2.4 


which has Serimene ee! 

h = A I.(u) » Bie (a) 
where I) and Ky are the modified Bessel functions. Since I (u)-»co 
and K (a) > 0 as u~moo, the solution to Eq. 2.4 which obeys the 
required boundary conditions is 


Q= Cr) Ko lu) 


where 
CS 


| 
—e rn 
C(a)=VCn) Fagier Ki (ELe a) 
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Under the Laplace transform, the acoustic pressure, radial component 


of particle velocity, and radial component of particle displacement 


become 


r so A vo (Vitam) 
ect tz) V (a) LotrWar  Y, (= attr ) 
Ue : a) ki (= J prt+ unt ) 

Z7(2) van (= Meier ) 


Ee | Via) Ki CE J ater ) 
£ (2) - oe Ky (> | peru” } 


Formal inversion yesults in expressions Like 
*¥ loo 


up = B® [We Ki(E fate | 
alte Ki (fate wr’) 


= tO 
where the contour is indented into the right half-plane around the 





s6 


fy 
A 
€ dap 


2 


singularities on the imaginary axis. Exact solutions to Eq. 2.5 


do not appear to be known, but certain approximations may be used 
to obtain approximate solutions for the leading portions of the 
propagating transient. 
: (15) 
For large values of argument, use of standard approximations 


for the modified Bessel functions results in the following asymptotic 


form 


Kn( =o) & ~V6 
wate [Ee Ft + Aum (Ss) 


ie 4 
+ Bram (e5)*+ Cam (<=) pk 


20 


1/ 


01 
— 33 6 OW Q CaN & 

B oreo = . eee — 

01° «12a * 728 7 Sa?) 
- ~~ 0. 

‘a> ~ g( 1-7) 

Ban = so tuBie DR cl 5. &\ 

- ie: i7@ re ~SBK PR : 


The symbol T represents the time of flight over the shortest 

path between the surface of the source and the field point. The 
restriction on s necessary to satisfy the approximation which led 
to Ea, 2,015 
OW 


|e oe 


qh 


Imaginary s greater than we in magnitude requires 


| 
- | = Wo 











where 


Wo = Wr 4+ (ees) 





is clearly a function of the source-radius and the cut-off frequency 


for the mode of interest. For a waveguide with pressure-release 
wa 
boundaries which is excited in the lowest mode, the quantity a 


a : : 
has the value Ws wich demonstrates the importance of the ratio 


of source-radius to layer-depth. As a result of the approximation 


which led to Eq. 2.6, Eq. 2.5 can be broken into two parts: 


t- Vo 
Ye om: = ( [2 [1 eA (Za)s Ba (se Jet Vda 


(z 
a 


+L U)ee 


+ \ V | tls) Pin. 
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Each portion contributes significantly to u- only when the exponent 
r 

of each is stationary. For large s (which corresponds to the 

leading portion of the signal), only the first integral is important 


and yields its major contribution when the exponent has a vanishing 


derivative ian 
So 


Pere 


d cass peels. : : 
S (at-Te)=O=t fpeteun® 


Thus, the requirement that ldo/fi | = UJasets an upper limit 
. mr 
ant | AwWw\e 
tas = | hs af ( =) 


on the time within which the first integral yields the major contribution. 





_ 1 
For a step input, for which V(s) a the radial component of the 


particle displacement becomes 
he | 2, ) 
Er 4 Ki ¢ CaN 
———ee = 
z (2) BD kK, (& (awe) 


Application of the asymptotic approximations of Eq. 2.6 leads to 


a a vs os 
Ht Ele — iu (de) or, 
iiveretome leads to 


© 
eet 


Ev a zl 


rw 
yy 
a, 


a 5 Len aan Sam(arg) 
wm! 
Leen 7 ) = Tama lara) 


+ 
ae 
i" e{r 
> 
\° < 
Ps 
; o 


where arg = We E> = Mes 


For a step input, the solution for the acoustic pressure becomes 


mel m\ ern | 
$< 7) sae 
ee Tay 2| VT. va) 1 Aoy ae a - 4] 


Sami larg) Bot (2) 2m (3) “Taalard? 





Oe, 


The quantity t' = t - T appearing in Eqs. 2.7 and 2.8 represents 
the time which has elapsed since the propagating transient signal 
raveling with velocity c first arrived at the position r. This 

quantity will be denoted as the "delayed time" 
For the leading portion of the propagating transient, one 
can apply a standard trigonometric asymptotic ee ee 
Tyg) = Vey cos by % ~ LE) . 
to the functions in Eqs. 2.7 and 2.8. Application of this approx- 
imation to the sums of Eqs. 2.7 and 2.8 leads to expressions of 


the types 


£ 
Se wil i GIT ~~ ( LT Q 
Tn a wsly 4 u 7 If lte™ 4 z 


f 
y 4nth oT & . 
fre conly Fo) lity Soa cosly Eas)?" 


Application of Eqs. 2.9 and 2.10 to the solution for particle 


av” <A casts it a ae form 


a +E zi YN + [EEE] §( 1-30) | 


IL ' ; Pl 
cos (2T1N "4 -T) san (Fc (2TTN - he | ; 
aw 
where 2TTN = y and No = a0 T — : 


The quantity N is the number of cycles of the transient generated 
in time t, and N. corresponds to the cycle where the approximations 
(and thus the theory) fail to apply. 

By combining the sine and cosine terms of Eq. 2.11, we find 
the expression for the radial component of particle displacement 


to be 


20 


Ty (teow 
Tt | — Le ca) sn (20N -& a + +) 


tz) mes 
where C(t) = (4-Run(S e\ + fan [* and 
tan 4 = = An dt Now 


1- Bu (NINa)* 


Application of Eqs. 2.9 and 2.10 to the solution for acoustic 


pressure casts it in the form 


le 2(4 7 i} iw Tr AL) cos (21N-& +) 23 


{C4-Bor(N/ne)? 12 + LAcItV/NS]? anc 


Aoi (N/Nad 
4 — Bor (N/Nwy © 


(5,14) 
Previous work has shown that for the same input as above 


where A(t) 


tan Y 


the exact solution for propagation of a transient in a duct becomes 
the leading term in the solution developed for the isovelocity layer. 
It should be noted that these duct terms can be retrieved from the 
trigonometric approximations by letting the amplitude factors go 

to unity and the phase angles go to zero. 

The solutions, using Bessel functions or asymptotic approxi- 
mations to the Bessel functions, for the transient response of an 
isovelocity layer have been developed with the restriction 

oe = Ga. Or Nye 

Computer programs have been written to evaluate the predicted 
waveforms for acoustic pressure and particle displacement, using 
either the Bessel function or trigonometric approximations (see 
Appendices A and B, respectively). These programs give printed 


and graphical outputs for the duct and layer terms. In addition, 


Zt 


these programs predict the axis-crossings and amplitude-modulation 
in the layer. Figures 2.1 and 2.2 illustrate the acoustic pressure 
and particle displacement, respectively, predicted by the two approxi- 
mations. Tables 2.1 and 2.2 give the calculated values of the predicted 
axis-crossings and amplitude-modulation for acoustic pressure and 
particle displacement. It should be noted that the first portion 
of the first half- cycle of the waveform predicted by the trig 
computer program is not accurate. To avoid a mathematical ambiguity 
at the leading portion of the waveform, its value has been set equal 
to zero. For comparison purposes, the predicted amplitude peaks 
were normalized to the value of the highest peak in the waveform. 
As can be Seen from a study of Tables 2.1 and 2.2, the difference 
between the two predictions is on the order of one per cent which 
is smaller than the experimental accuracy encountered. Therefore, 
the predictions from the computer program using trigonometric 
approximations will be used for all comparisons with experimental 
data. 

The approximate solutions for the prepagation of a transient 
in an isovelocity layer have been developed, and computer programs 
to evaluate these solutions and predict the axis-crossings and 
amplitude-modulation in the layer have been written. 

Transfer Characteristics of Mylar Transducers 

To have some measure of confidence in the received waveforn, 
it was decided to investigate the transfer characteristics, under 
steady-state conditions, for Mylar transducers. This investigation 
will cover the radiation impedance at the surface of the transmitter, 


r =a, and the radiation impedance of a receiver. 


Pe 


wo 
@ 
© 


© | 
: 
EF G0 O08 \/ 6 oe we t' 
5 
aa 
Be) 
a 
a8 
| Computed using Bessel functions 
: 
B / 
B EP OE c ~, ae 
iB 
oe \ : : 
t 
all \ 
aT Jomputed using trigonometric approximations to 


Bessel functions 


Acoustic Pressure p vs. Delayed Time t' as Predicted by Computer 
Programs for the Conditions r=91.) cm, a=.00 em and d=8.00 cm. 


Figure 2.1 
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Table 2.1 


AXIS CROSSINGS AXIS CROSSINGS AMPLITUDES OF AMPLITUDES OF 
USING BES USING TRIG MAXIMA, BES MAXIMA, TRIG 
TIME IN USEC TIME IN USEC NORMALIZED NORMALIZED 
1.000 1.000 
1.45 139 
0.632 0.647 
oo 7.46 
0.470 0.478 
le .32 L8e26 
Cro20 Owe96 
339098 33), Sill 
0.389 0.344 
pi2e99 S259 1. 
0.303 0.308 
76.17 76.08 
OaZ 75 0.280 
LOZee2 102.64 
0.253 0.258 
132633 132421 
07255 05259 
164.57 164 .42 
0.212 On224 
Tors 1982°96 
0.205 Os 10. 
Zor 2335¢o04. 
0,193 0.198 
274,07 Zio ee 
OnhkBZ 0.186 
Sl og0 313.65 
0.171 0.8/7 
Soe 354.81 


Comparison of computed axis-crossings and amplitudes of maxima for 
acoustic pressure using Bessel functions and trigonometric approxi- — 
mations to Bessel functions for the conditions r = 91.4 cm, a = 4.00 cm, 


and d = 8.00 cm. 
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Ce) 


a i a 
DO} o, f\ oy, x aa 4! 


es » arbitrary units 
096 


7) Computed using Bessel functions 


© 


B04 


G3s 





Computed using trigonometric approximations to 
Bessel functions 


= 4 arbitrary units 


-003 


Radial Component of the Particle Displacement e vs. Delayed Time t' 
Predicted by Computer Programs for the Conditions r=91.8 cm, a=6.00 cm, 


and d=8 .00 CMe 


Figure 2.2 
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Table 2.2 


AXIS CROSSINGS AXIS CROSSINGS AMPLITUDES OF AMPLITUDES OF 


USING BES USING TRIG MAXIMA, BES MAXIMA, TRIG 
TIME IN USEC TIME IN USEC NORMALIZED NORMALIZED 
0.474 0.475 
Sa 3292 
0.701 0 698 
12.48 12766 
0.844 0.841 
25.94 25moo 
0.95835 0,2 
43.82 43.79 
Over 0.981 
Goa 5 65.64 
1.000 1.000 
Wes ola 
0.9) Ong26 
120.16 Lhome? 
0.974 Ons75 
Le. Oo IStor 
0.941 0.943 
186.05 185.65 
O2901 Omvu>5 
222.44 22062 
0.857 0.860 
ZO0%7 3 260.24 
0.811 0.814 
300.68 300% 
0.764 0.768 
sae 07 341.50 


Comparison of computed axis-crossings and amplitudes of maxima 
for the acoustic particle displacement using Bessel functions and 
trigonometric approximations to Bessel functions for the conditions 


r = 91.8 cm, a = 6.00 cm, and d = 8.00 cm. 
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In a steady-state condition, the acoustic potential for a 
cylindrical source in a waveguide can be expressed as 


P, = (2) Rlr) ots 2.14 


The time independent wave equation becomes 


4° | ddr 
eet, a + ide =O. 


The separation of variables suggested in Eq. 2.14 leads to the 


/ W 
value of the separation constant, ke +t oe where the phase velocity 


P 
a is defined as 
a 
Ch = Pies 2.15 
gerbe , 
t 4- (#8) 


Thus, the wave equation for the radial component of the acoustic 


potential becomes 


d*k Pea. Ee 


“ie 2 yr ar 


with the general solution 
R(r)= A Xo (ker) s 6 Ne (uy), 2.16 


For an outward traveling wave, Eq. 2.16 becomes 


(2) 
R (r)= Al J. (ker) -b To (ky r)| = Alo (&rr). 2.17 


The acoustic potential, acoustic pressure, and the radial component 


of the particle velocity bec ug 


dr yo “sr 
ele) ~ (% tee 


qa) 


os ve A oe) (se ‘ 


(é 


r Ciena 
-iwe 0 ee )e 
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Using Eq. 2.15, the specific acoustic impedance Z becomes 


z= fo = dec, He tom as 


PIE: 


where arg = = \ \- (~2)° 


Then, the radiation impedance Z_ is given by 


Z eee) bio’ ® sd psn 
Z-=AL - re Wie) Caray 


where A is the total radiation area of the transducer. Thus, the 


radiation impedance at the surface of the transducer, r =a, is 


given by 
, Cd (wer [4 Wn \2 
Z L = a \=& sabe Lal 
y -_ 


CB yt (STB) © 


The behavior of Eq. 2.19 as a function of frequency will be 
investigated. An estimate of the maximum frequency encountered 
in a shallow-water layer can be derived by recalling that the 
highest frequencies in the transient occur in the first half-cycle. 
The following relation concerning the time for the maximum frequency 


can be derived: 


Wn Vt man 7 Te =2TIN=Tl or 


= 
Cte = eS 4 ee where 2a V° 


The relation for the instantaneous frequency 


Wt want ee 
oe bee 


can then be manipulated into the desired result 


Winn = i + (ey ° 
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For r-a = 100 cm and d = 8 cn, 


Ww 650 khz. 


max 
The cut-off frequency for the first mode and the above approximations 
1s 
w. = 59 khz. 
n 
Thus, the range of frequencies expected to be encountered are 


wi aw <w The expected frequency must always be greater than 


max- 
the lower limit, for the assumption of a steady state condition no 
longer holds at w = wo: 

Since the Hankel functions are complex numbers, they can be 
expressed as a modulus and a phase angle. Thus, Eq. 2.19 can be 
rewritten as 

BCs e GS Mod ao 

A f 
where Mod is the modulus and 9 is the phase angle. Figure 2.3 
shows the dependence of the phase angle and modulus on the argument 
of the Hankel functions. For small values of the arguments, the 
radiation impedance is strongly dependent on the frequency. How- 
ever, when the value of the argument is about Tl which corresponds 
to a frequency within 0.5 per cent of the cut-off frequency, the 
radiation impedance can be approximated by 2@c Thus, for the 
frequency range of interest, the radiation impedance is frequency 
independent. 
The open-circuit mechanical impedance of a Mylar eeateducer = 
is given by Z =r + i(wm - s/w) where r is the transducer mechanical 
resistance and s is the transducer mechanical stiffness. For the 


frequency range of interest, the mechanical resistance is considered 


to be negligible in comparison with the reactive part of aa° Under 
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(16) 


the assumption that the mechanical stiffness is given by 
c= PAS, where i. is the atmospheric pressure and f is the 
effective depth of the air cushion behind the Mylar diaphragm, the 


equation for Z can be rewritten as 


m 
f= CA we t - Rae 
Lm LA L Cavin 4 
where t is the thickness of the Mylar. With the nominal values 
= 1.395 Silene 
Pry Lax ‘ 


ey Semon cm, 


Pa Polo x toe Neem’, and 


Q= 2.00 x 10> em, 


the values of Zs for the frequency range of interest can be 


calculated, Figure 2.4 shows the dependence of the magnitude of 


a 
rk Se frequency. 
2 32 (16) s 
The velocity of the diaphragm of a Mylar transducer is 
given by 
Rov Co 


ee weC. (Zwrtr) -Or+iwlal@merer) 


where Co is the blocked diaphragm electrical capacitance, v is 
the signal voltage across the transducer, § is the turns ratio, 
and R is the internal resistance of the power amplifier. The 


turns ratio @ is defined as 


_ CoV 
PS. ee 


where ve is the external polarizing voltage and XG is the position 
of the Mylar film when the polarizing voltage is applied. A 


comparison of the magnitudes of the terms in the denominator of 
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is proportional to the particle displacement. These conclusions 
only hold for the leading portions of the received transient where 
Z. can be replaced by ee 

The effect of the transfer characteristics of Mylar transducers 
has been found to be dependent on the frequency and the radius of 
the transducer. Given a transient signal which is the mirror of 
the input voltage, the following is a summary of the effects that 
would be observed as a result of the radiation impedance: 

(a) since the radii of the family of transmitters were large, 
the response is essentially flat throughout the region of interest; 
thus, no effects should be observed. 

(b) because of the small radius of the receiver, the effect 
becomes more pronounced than for the transmitters. One should 
observe a reduction in the size of the peaks of the leading portion 
of the transient and should observe an enhancement of the peaks as 
the local apparent frequency approaches the cut-off frequency. A 


similar phase angle effect should be observed. 
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3, EXPERIMENTAL APPARATUS 


The experimental apparatus in this research project consisted 
of a tank, a family of cylindrical Mylar transducers, a barium 
titanate array receiver, and the associated electronics equipment. 
The laboratory setup is shown in Fig. 3.1. 

Tank 

The water layer was contained in a tank constructed of 5/8 
inch marine plywood. The dimensions of the tank were 6 ft x 6 ft x 
1 ft. The plywood was covered with Varathene to provide water- 
proofing and lined with 1/4 inch polyethylene foam which was assumed 
to provide the desired pressure-release bottom boundary of the 
layer. The polyethylene, not being waterproof, was coated with 
liquid Neoprene paint to prevent deterioration of the pressure- 
release surface. The bottom of the water layer was maintained 
level by means of adjustable feet installed in the legs of the 
Supporting platform. Since the walls of the tank were also Lined 
with polyethylene, the tank could be excited as a pressure-release 
cavity, which provided one method of determining the depth of the 
layer. 

Transducers 


Previous on 


had indicated that Mylar dielectric 
transducers were well-suited for use in layers. 

A family of eight Mylar transducers and one barium titanate 
linear array were constructed for use in this investigation. The 


radii of the Mylar transducers were 10 cm, 6 cm, 4 cm, 2 cm, 1 cn, 


0.5 cm, and two at 0.25 cm. Each had an active radiating face 


oF 


8 cm in height which was the maximum depth of the layer. These 
various transducer dimensions were chosen to provide a large range 
of values for the ratio N/N,- The bérium titanate arsay was 8 cm 
high and 1/8 inch in diameter. 

The Mylar transducers were constructed using techniques 


(16) 


described in , modified as follows: 

(a) The backplates of the larger transducers were constructed 
from 4% inch thick aluminum cylinders. The backplates for the 
smallest transducers were fabricated from % inch diameter copper 
tubing. 

(b) The cylinders were capped at the top and bottom by 
lucite plates. 

(c) The Mylar film was cut to the desired dimensions and 
attached to the aluminum cylinders with thin strips of double-stick 
tape, one piece around the top, one around the bottom, and one 
piece vertically from bottom to top. 

(d) The inner electrical connection was made from the center 
pin of a BNC connector mountzd in the center of the top lucite 
plate to the inside of the backplate. The outer conducting path 
from the metallized side of the Mylar film to the BNC was completed 
by covering the top lucite plug with aluminum foil which was lLapped 
over the Mylar film. To ensure good electrical continuity, silver 
paint was applied to the aluminum-Mylar interface. To decrease the 
possibility of electrical pickup, the bottom lucite plate was also 
covered with aluminum foil. 

(e) When the silver paint had dried, the entire transducer 
was covered with liquid neoprene to provide the necessary water- 


proofing. 
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The pressure-sensitive barium titanate was constructed using 
eight 1/8 inch by 1/8 inch cylinders mounted on several strands 
of twisted wire, which provided the inner electrical connection. 

To provide the damping necessary to reduce the effect of Longitudinal 
resonances, the spaces between the cylinders were filled with strips 
of rubber. The outer electrical connection was provided by a wire 
which was coiled around each cylinder. This probe proved to be 
somewhat insensitive and was used only for comparison purposes. 

When the electrical load on the Mylar receiver (see Fig. 2.5) 
was terminated by the 1 megohm impedance of the voltage amplifier, 
the receiver was sensitive to particle displacement. Therefore, 
the Mylar receivers were used to measure the particle displacement, 
and the barium titanate array was used to meaSure the acoustic 
pressure. Transmitters and receivers are shown in Fig. 3.2. 

Electronic Equipment 

Figure 3.3 shows the block diagram of the equipment; Table 
3.1 identifies the components in the figure. The unit-pulse generator 
was used for the step input. 

The received signal was sent from the receiver through the load 
and bias network, amplifier and filters to trace "B" of the dual- 
beam oscilloscope. The output of the power amplifier could be 
displayed on trace "A", Measurements of the received waveform axis 
crossings and amplitudes were made directly from the oscilloscope 
display, or photographed for later analysis. The electronics were 
triggered by the time-mark generator which also triggered the time- 
base of the oscilloscope. The intensity of the ''B'' trace was 


modulated by a 10 usec marker available from the time-mark generator. 
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VEG 


PWR AMP 


TMG 


AMP 1 


SKL FIL 


RC FILL 


RC FIL2 


SCOPE 


LBN 


Table 3.1 


Electronics Equipment 


Unit pulse generator, General Radio type 1217-B 

Power amplifier, Hewlett Packard 467-A 

Time Mark Generator, Tektronic type 180A 

Voltage amplifier, Tektronic type 1121A 

Electronic band pass filter, SKL model 202-D 

Homemade variable low-pass RC filter with 400 khz roll-over 
Homemade high-pass RC filter with 2 khz roll-over 

Dual beam oscilloscope, Fairchild model 777 


Load and bias netwrk shown in Fig. 3.4. 
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Source and receiver were carefully positioned so that the distance 
between adjacent sides, r-a, was always the same. This insured 
that the trace modulation would be properly positioned with 
respect to ye inception of the received transient signal. The 
load and bias network shown in Fig. 3.4 served to supply the D.C. 
voltage bias needed by the Mylar transducers without allowing 
this voltage to appear elsewhere in the circuit and to provide 
the proper electrical load for the receiver. The 22 ohm resistor 
in the transmitter circuit was needed to Suppress parasitic 
oscillations caused by the transient excitation of the capacitance 
of the transmitter and the residual inductance of the power 
amplifier. 

The rise time of the system was experimentally determined 


to be about 0.5 usec. 
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4, EXPERIMENTAL INVESTIGATION 


The purposes of the experimental investigation were to (a) observe 
and analyze the received acoustic pressure and radial particle displace- 
ment waveforms in the layer, (b) compare these waveforms with the 
theoretical predicted waveforms produced by the computer programs, and 
(c) vary the value of the ratio N/N. over a large range to investigate 
the validity of the assumed asymptotic approximations previously dis- 
cussed, Four basic parameters characterized each experiment: the 
range r, the transmitter radius a, the layer depth d, and the speed 
of sound, c. 

Experimental Procedures 

Before each run,it was decided which parameter would be varied. 
Generally, the time of flight T and the depth of the layer d were 
kept constant while the source radius a was varied throughout its 
range. The 0.25 cm radius Mylar transducer was used to record the 
received particle displacement, and the barium titanate linear array 
was used for the acoustic pressure receiver. The temperature of the 
water was measured to determine the speed of sound in the eee 
The depth was determined by averaging several physical measurements 
taken over the ensonified area. Another method of determining the 
depth was to excite the tank as a pressure-release cavity. The 
cut-off frequency could then be determined by observing the received 
Signal on an oscilloscope and varying the frequency until the signal 
disappeared, Once the cut-off frequency was known, the depth could 
be calculated, using the formula wae a . Because of the irregularities 


observed in the depth, this method proved to be no more accurate than the 
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actual measurement uSing a ruler and was consequently abandoned. 
The intensity of the trace of the oscilloscope was modulated by 

a 10 usec marker available from the time-nark generator. The 
source and receiver were carefully positioned so that the distance 
between adjacent sides, r-a, was always the same and such that 

the reception of the received transient and one of the time marks 
always coincided. This insured that when changing transmitters, 
identical positioning and time of flight could be achieved. 

The received waveforms were photographed from the oscil- 
loscope display. For the step input, several photographs were 
taken of each waveform. An expanded time scale was used when 
photographing the leading portion of the waveform which conteined 
the high frequencies. For a typical experimental run, three 
photographs of each waveform were made utilizing, for example, 

20 usec/cm time scale. By varying the delay time vernier, the 
entire waveform could be displayed and photographed for Later 
analysis. By careful measurement of the position of each exis 
crossing with the neighboring time marks, the relative phase of 
each axis crossing could be compared between the photographs taken 
for different source radii. The amplitudes of the peaks of the 
waveforms were read with the help of dividers and a scale. Sample 
photographs are shown in Fig. 4.1. It should be noted that the 
anomalous behavior of the received waveforms 4s evidenced by 

Fig. 4.1 departs drastically from the smoothly-varying computer 
predictions of Figs. 2.1 and 2.2. This anomaly will be discussed 


in more detail later. 
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bserved Waveforms 


For each experimental run, a corresponding computer run was 
made uSing the same values for r, a, c, and d. The computer 
programs were written so that their cutput for e2ch run included 
graphs of the predicted waveforms, times of axis crossings for 
the layer solution te and for the duct solution tom? and amplitudes 
of the peaks and the time of occurrence of these peaks for the 
layer solution. In this manner, the theoretical values of the 
times of axis crossings and of the amplitude peaks could be 
compared with those values experimentally determine. 

The time of thé axis crossings were compared in the following 


(5) 


manner: Previous work suggested that the time of 2ach axis 
crossing should be the same as the theoretically predicted t 
within the validity of the epproximetions used. As w2s shown in 
section 2, the theory was valid only as long as N/N. l. There- 
fore, the convenient parameters for analyses were the time of the 
axis crossings and the ratio N/N.. In ord@r to reduce the size of 
the graphs, the time coordinate was reduced by the amount S predicted 
by the computer.. The time coordinate signified by At! was pletted 
against the ratio N/N,. Figures 4.2 to 4.4 are representative 
examples of these graphs. The solid line represents the values of 
Ces 5 ty and the circled dots represent the values of 2 - ta 
Were the theory to have been completely correct and the experiment 
accurate enough, the data would have fallen on the X-exis (which 
represents t - t)- As can be seen from Figs. 4.2 to 4.4, the 


data were not in agreement with the theory. The effect of the 


phase angle demonstrated in Fig. 2.6 was taken into account and 
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plotted on Figs. 4.3 and 4.4. As can be seen, this did not 
improve the agreement between the theory and the data. 

The peak-amplituce data were also analyzed. The experimental 
amplitude peaks were first normalized to the value of the highest 
peak in the waveform. The predicted peeks were normalized in the 
course of the computer programs. The ratio of these normalized 
quantities was then piotted as a function of the delayed time of 
the predicted amplitude peaks for each run. Exact agreement in 
the shape of the transient would result in a line of zero slope 
passing through the point (1,0). Tne amplitude effect discussed 
in Section 2 for a Mylar transducer used as a receiver was also 
plotted. Figures 4.5 to 4.7 show plots tor various values of r, 
a, and d. Once again, this effect did not improve the agreement 
between the theory and the data. 

Error Analysis 

The rise time of the experimental system was determined to be 

about 0.5 usec. Therefore, instead of being a perfect step-function 


input, the velocity at the source was initially described by 


V(t) : 1 > a “ae 


where of = ILS 10° oi is the reciprocal of the rise time. 


(3) 


Previous work has indicated that the rise time of the system 
should introduce distortion into the recéived signal fcr at most 
3.0 usec of delayed time t'. Since this time is less than that 


for the first maxima, the effect of the rise time on the received 


Signal can be ignored in this investigation. 
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It has been shown that the velocity of the Mylar face decays 
to 88 per cent of its maximum value after a time of 360 gee? 
This effect is not large enough to eccount for the tailing off 
of the amplitude peaks of the received waveforms seen in Figs. 

G25 tora. 74 

Uncertainties in r, a, c, and d were the possible sources 
of error in this investigation. The effective radius of each 
transmitter was measured to the nearest 0.5 mm, which could have 
introduced an uncertainty of + 1.0 per c#nt in the value of a 
for the 6.0 cm radius source and + 5.0 per cent in the value of a 
for the 1.0 cm radius source. Measurement of the temperature of 
the water to the nearest half-degree restricted the systematic 
error inc to within + 0.5 per cent. Determination of the time- 
of-flight T of the received signal form the source to the receiver 
produced a possible systematic error of + 0.5 per cent, which in 
turn resulted) m an erronan r of G 0.>5aperecen-. 

The largest experimental error was encountered in determining 
the depth of the layer. Observed irregularities of + 1.0 mm in 
the bottom surface of the tank could introduce an error in the 
local effective depth of the layer. Since the averaged layer depth 
used as an input to the computer programs was only accurate to + 1.0 
mm, the systematic error introduced in w, was on the order of + sO 
per cent for a depth of 8.0 cm. 

The error in determining the times of the axis crossings fron 
the data was effected by the accuracy of the oscilloscope time-base 


which was calibrated to within 0.3 per cent. With the trace of 
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the oscilloscope modulated with the 10 usec marker as an additional 
time reference, this error was reduced even more. Therefore,a 
random error of + 1.0 per cént could be expected when one takes 
into account the accuracy of reading the oscilloscope scale. 

The errors in the computed values of t were influenced 
predominately by the error in w_, since it appears in the argument 


2 
of the Bessei functions (we fees 1 ) from which the t are 


oS 


computed. It has been shown that the argument of the Bessel 


functions could be rewritten as 


Wi Fen (tint ZT) nai ae 


where Ce are a set of constants which yield zero values for the 
amplitude of the particle-displacement solution. For the leading 
portion of the signal, we have 

wntel nt ~ Cun 


and es. 
caer nna eeeeremeee en 
Con Meee 
Therefore, the + 3.0 per cént error in we could lead to a + 6.0 
per cent error in the computed values of t 
To reduce the scatter of the data, a one-parameter fit was 
tried. It was decided to determine an effective cut-off frequency 
and use this es an input to the computer programs instead of the 
layer depth. This was accomplished by determining the cut-off 
frequency for the values of ae and the arguments of the trig 
function approximaticns developed in Section 2. At an axis 


crossing, the argument was equal to an integer multiple of TT 


wnt alte tel) - % +Ylewmll, 


in 
ae 


A computer program was written to evaluate the cut-off frequency 
for each axis crossing. Then, these cut-off frequencies were 


' to be used in 


averaged to give an "effective cut-off frequency' 
the computer prcgrams to recompute the time for each axis crossing. 
Figure 4.8 shows the effect of this one-parameter fit. As can be 

seen, it did not improve the scattering significantly and was not 

used. 

Error flags on Figs. 4.2 to 4.8 indicate the relative magnitudes 
of errors in the experimental and computed values of the axis crossings 
and peak amplitudes. However, these error flags are not large enough 
to explain the scatter of the data or the anomaly which was present 
(see Fig. 4.1). The theory could not account for this anomaly so 
various experiments were conducted to investigate its character and 
source of origin. 

Firstly, it was suggested that irregularities in the bottom of 
the tank might be a possible cause of the anomaly. Therefore, the 
relative positions of the transmitter and the receiver in the tank 
were changed. The values of T and d were always kept constant. Four 
contrasting positions were chosen; one near the left hand side of 
the tank, one in the center of the tank, one adjacent to a seam in 
the bottom near the right hand side, and one where the transmitter 
and receiver straddled the seam. As can be seen in Fig. 4.9, within 
the experimental accuracy of the investigation, positioning on the 
bottom had no effect on the received waveferm. 

Secondly, it was decided to monitor the current drawn by the 
transmitter when loaded by the water medium and compare this with 


the current predicted assuming a loading. The current was 
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experimentally determined by meaSuring the current in the two legs 
o£ the circuért~shownin~ Fig +.t0e7-—The=current present in each 

leg was amplified by a differential amplifier and then the difference 
displayed on the oscilloscope. Photographs of the received signal 
could then be analyzed to determine the current drawn as a result 

of acoustic radiation into the medium. By analyzing the equivalent 
circuit of a Mylar transducer used as a CAGE Gaines as shown 

in Fig. 4.11, the predicted value of the current drawn as the result 
Oot the aS ee loading can be determined. For this assumed 
loading, the Serie drawn by the mechanical leg of the source 
should be flat. Typical values for the 6 cm radius source were 

near 6 uamp. The measured values of the current had rapid oscil- 
lations of nearly 1 uamp peak-to-peak amplitude, which indicated a 
significant fluctuation from the flat response predicted. This 
casts doubts on the validity of the assumption of me loading used 
in the developments of Section 2. 

Lastly, it was suggested that the polyethylene foam was not a 
perfect pressure-release surface and that transmission through it 
into the wood bottom was the cause of the anomaly. In order to 
investigate this possibility, two additional layers of polyethylene 
were glued to the bottom of half the tank. Then runs were made on 
the single-layered bottom and then on the three-layered bottom. 

The quantities r, a, and d were kept constant for the corresponding 
runs. Figures 4.12 to 4.14 indicate the predicted waveform, the 
experimentally determined waveform on the three-layered bottom, and 
the experimentally determined waveform on the single-layered bottom 


for three values of a: 10, 6, and 4 cm. The vertical and horizontal 
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22 ohms 0.5 uf | Transducer | 







Pwr Amp 
Output 


C. in leg B was matched to value of Co of transducer (leg A). 


Circuit used to Monitor Current Drawn by a Mylar Transducer as a 
Result of Acoustic Radiation into a Water Mediun. 


Figure 4.10 
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Equivalent Circuit of a Mylar Transducer used as a Transmitter. 


Figure .11 
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(b) Received particle displacement waveform on three-layered 
bottom 
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(c) Received particle displacement waveform on single-layered 
bottom 


Bottom Layer Investigation for Three-Layered Bottom and Single-Layered 
Bottom using 10.00 cm Source. 


Figure 4.12 
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(c) Received particle displacement waveform on single-layered 
bottom 


Bottom Layer Investigation for Three-Layered Bottom and Single-Layered 
Bottom using 6.00 cm Source. 


Figure ).13 
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Bottom Layer Investigation for Three-Layered Bottom and Single-Layered 
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Figure 1h 
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scale settings on the oscilloscope were kept the same for each 
set of runs. As can be seen, there is some evidence of a bottom 
effect, but there is no clear correlation between the three runs. 
However ,this evidence is sufficient to doubt the credibility of 
the assumption of the perfect pressure-release characteristics of 


the 0.25 inch polyethylene foam. 
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5. RESULTS AND CONCLUSIONS 


The investigation has shown that for N/N_S 1 the predicted 
and experimentally observed shifts in axis crossings should agree 
to within + 6 per cent if one takes into account just the known 
systematic errors. The correlation between the predicted and 
observed amplitudes should be of the same accuracy. However, 
the difference between the theory and data is greater than that 
expected, 

The theory as developed in Section 2 could not account for 
the anomaly which was present in all runs, but the qualitative 
agreement between the observed and predicted waveforms and the 
general agreement with the solutions developed in other pea ee 
lead one to believe that the theory presented in Section 2 is correct. 

It is concluded that the experimental setup used is not accurate 
enough to give results needed to confirm the theory. 

It is recommended that better electronic equipment and transducers 
should be investigated. Further studies should be conducted to check 
the validity of the assumed pressSure-release characteristics of the 
polyethylene foam, More investigations are needed to try to determine 


the cause of the observed anomaly. 
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